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^ ! Abstract 

I A theoretical study of some electrodynamic features of a region close to a slowly-rotating mag- 

• netized relativistic star is performed. To be a little more specific, based on the solution-generating 

. method given by Wald, the magnetic fields around both uncharged and (slightly) charged rela- 

^ ■ 

tivistic stars have been obtained. Particularly for a charged relativistic star, again following the 

o ■ 

^ I argument by Wald, the star was shown to gradually accrete charge until it reached an equilibrium 

■ value Q = 2BqJ. This value of the equilibrium charge seems to be generic as a rotating black hole 



is known to accrete exactly the same amount. Although these results are equally relevant to all 



H ' species of slowly-rotating relativistic stars, we particularly have the rotating neutron star in mind. 

As such, it would be of some interest to attempt to make contact with a real pulsar case. Thus, 
we discuss how many of the theoretical results obtained in the present work can be carried over to 
a realistic, general relativistic description of a pulsar's magnetosphere. 
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I. INTRODUCTION 



In the present work, we perform a theoretical study of some electrodynamic features of a 
region close to a slowly-rotating relativistic star. Here, by relativistic star, we mean a highly 
self-gravitating compact object. Particularly, we have a rotating neutron star in mind. Thus, 
we first begin by providing the rationale for treating the vicinity of a spinning compact 
object, such as a rotating neutron star, relativistically as a non-trivial curved spacetime. 
Since the pioneering proposals by Gold and by Pacini j^, it has by now been widely 
accepted that, indeed, pulsars might be rotating magnetized neutron stars. Nevertheless, 
since then, nearly all studies on pulsar electrodynamics have been performed by simply 
treating the region surrounding the rotating neutron stars as flat. This simplification may 
be sufficient just to gain a rough understanding of the origin of the pulsars' radiation. 
However, nearly thirty years of studying pulsar electrodynamics have passed, and now, it 
seems we should look into the situation in a more careful and rigorous manner. 

In the present work, we take one step forward in this direction. We particularly aim 
at determining the amount of charge on and the structure of the magnetic field around a 
slowly-rotating neutron star by treating its exterior region, literally, as a curved spacetime 
and employing the algorithm of Wald [2| . Indeed, thus far, only a handful of earlier studies 
concerning the electrodynamics of rotating neutron stars have appeared in the literature, 
and they include the works by Muslimov et al. 4] and by Prasanna and Gupta 4]. Although 
their motivatioUjapproach, and tools are completely different from those of our present work, 
RezzoUa et al. 'q] also considered analytic solutions of Maxwell equations in the internal 
and external regions of a slowly-rotating magnetized neutron star. 

Indeed, we now have a considerable amount of observed data for various species from radio 
pulsars to (anomalous) X-ray pulsars Even if we take the oldest-known radio pulsars 
for example, it is not hard to realize that these objects are compact enough that they should 
really be treated in a general relativistic manner. To be more concrete, note that the values 
of the parameters characterizing typical radio pulsars are ro(radius) ~ lO^(cm), M(mass) ~ 
I.4M0 ~ 2x 10^^[g), r(pulsation period) ~ 10"'^ — l(sec), and i?(magnetic field strength) ~ 
10^^ (G). Thus, the Schwarzschild radius (gravitational radius) of a typical radio pulsar is 
estimated to be rsch = 2GMq/c^ ~ 3 x lO^(cm) (where G and c denote Newton's constant 
and the speed of light respectively) ; hence, one ends up with the ratio ro/rsch ~ 10^(cm)/3 x 
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lO^(cm) ~ 3. This simple argument indicates that, indeed, even the radio pulsar (which 
is perhaps the least energetic among all of pulsars) is a highly self-gravitating compact 
object that needs to be treated relativistically. Of course, a "relativistic" treatment here 
means that the region surrounding the pulsar, namely, the pulsar magnetosphere, has to be 
described by using a curved spacetime rather than simply a flat one. Then, the question 
now boils down to this ; what would be the relevant metric for describing the vicinity of a 
rotating neutron star ? Although it does not seem to be well-known, fortunately, we have 
such a metric for the region exterior to slowly-rotating relativistic stars, such as neutron 
stars, white dwarfs, and supermassive stars, and it is the one constructed long ago by Hartle 
and Thorne j^. Thus, the Hartle-Thorne metric is a stationary axisymmetric solution to the 
vacuum Einstein equation and in the present work, we shall take this Hartle-Thorne metric 
as a relavant one to represent the spacetime exterior to slowly-rotating neutron stars. 

As we mentioned above, here in the present work, we are particularly interested in the 
exact theoretical determination of the amount of charge on and the structure of the magnetic 
field around a slowly-rotating neutron star by treating its exterior region literally as a curved 
spacetime. As such, it seems relevant to introduce a related recent work which, as a matter 
of fact, has been one of the motivations for this study. In our earlier work P], we explored 
how much of the magnetic flux could actually penetrate the horizon of a rotating black hole 
in association with the Blandford-Znajek mechanism for rotational energy^xtraction 
from a Kerr hole. In doing so, we made use of the algorithm suggested by Wald 3] to figure 
out the amount of charge on and the structure of the magnetic field around a rotating black 
hole. The results presented there for the case of a rotating black hole can, therefore, be 
compared with those that will be presented in this work for the case of a rotating neutron 
star because, as stated, we shall in this work employ the algorithm of Wald to determine 
the amount of charge accretion on and the structure of the magnetic field around a rotating 
neutron star, as well. Thus, it should be interesting to read the present work in parallel 
with our earlier work on Kerr black hole electrodynamics. 

Speaking of the Blandford-Znajek mechanism, perhaps it would be relevant to briefly 
sketch their idea here as it has a significant impact on the physics underlying both the black 
hole and the (rotating) neutron star's magnetosphere. According to the Blandford-Znajek 
mechanism the way some of the hole's rotational energy is carried off by the magnetic 
field lines threading the hole can be summarized as follows : Angular momentum and ro- 
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tational energy transport is achieved in terms of the conservation of the "electromagnetic" 
angular momentum and the energy flux flowing from the hole to its nearby fleld particularly 
when the force-free condition is satisfled in the surrounding magnetosphere. It is also amus- 
ing to note that the Blandford-Znajek mechanism can be translated into a simple-minded 
circuit analysis. In this alternative picture, the current flowing from near the pole toward 
the equator of the rotating black hole geometry and the poloidal magnetic fleld together 
generate a "magnetic braking" torque that is antiparallel to the hole's spin direction. As a 
result, the rotating hole, which is a part of the circuit, spins down, and its angular momen- 
tum and the rotational energy are extracted. For more details, we refer the reader to the 
original work of Blandford and Znajek and some com pre hensive review articles jll| on 
this topic, including the "membrane paradigm" approach 1^ . 



II. THE HARTLE-THORNE METRIC FOR THE REGION EXTERIOR TO 
SLOWLY-ROTATING NEUTRON STARS 
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The Hartle-Thorne metric j8[ is given by 
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with M, J and Q being the mass, the angular momentum, and the mass quadrupole moment 
of the slowly rotating neutron star, respectively, P2(cos^^) = (3cos^^^ — l)/2 being the 
Legendre polynomial, and being the associated Legendre polynomial, namely. 
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As is well-known, the only known exact metric solution exterior to a rotating object is the 
Kerr metric . Thus, it would be worth clarifying the relation of the Hartle-Thorne metric 
br slowly-rotating relativistic stars given above to the Kerr metric. As Hartle and Thorne 
^, take the Kerr metric given in Boyer-Lindquist coordinate and expand it to second order 
in angular velocity (namely, the angular shift P'^) followed by a coordinate transformation 
in the fr, ^)-sector, 
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where a = J/M. Then, one can realize that the resulting expanded Kerr metric coincides 
with the particular case Q = J'^/M (with Q being the mass quadrupole moment of the 
rotating object) of the Hartle-Thorne metric. Therefore, in general, this Hartle-Thorne 
metric is not a slow-rotation limit of the Kerr metric. Rather, the slow-rotation limit of 
the Kerr metric is a special case of this more general Hartle-Thorne metric. As a result, 
the Hartle-Thorne metric with an arbitrary value of the mass quadrupole moment Q can 
generally describe a slowly-rotating neutron star of any shape (as long as it retains the 
axisymmetry) . 

Next, since we shall employ in the present work the Hartle-Thorne metric to represent the 
spacetime exterior to slowly-rotating neutron stars, we would like to carefully distinguish 
between the horizon radius of the Hartle-Thorne metric and the actual radius of the neutron 
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stars. Indeed, one of the obvious differences between tlie black liole case and tlie neutron star 
case is the fact that the black hole is characterized by its event horizon while the neutron 
star has a hard surface. Since this solid surface of a neutron star (which we shall henceforth 
denote by tq) lies outside of its gravitational radius, which amounts to the Killing horizon 
radius of the Hartle-Thorne metric, at which A = (1 — 2M/r + 2J^/r^) = in Eq. (3), 
we have tq > th- In the present work, however, we shall never speak of the Killing horizon 
of the Hartle-Thorne metric as it is an irrelevant quantity playing no physical role. 

Now, we turn to the choice of an orthonormal tetrad frame and we shall particularly 
choose the zero-angular-momentum-observer (ZAMO) ^| frame, which is a sort of fiducial 
observer (FIDO) frame. Generally speaking, in order to represent a given background ge- 
ometry, one needs to first choose the coordinate system in which the metric is to be given 
; next, in order to obtain the physical components of a tensor (such as the electric and the 
magnetic field values), one has to select a tetrad frame (in a given coordinate system) to 
which the tensor components are to be projected. As is well-known, the orthonormal tetrad 
is a set of four mutually orthogonal unit vectors at each point in a given spacetime. And 
these unit vectors give the directions of the four axes of a locally-Minkowskian coordinate 
system. Such an orthonormal tetrad associated with the Hartle-Thorne metric given above 
may be chosen as = e^dx^ = (e°, e^, e^, e'^). 
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A local, stationary observer at rest in this orthonormal tetrad frame has a worldline 
given by {dr = 0, dQ = 0, [d^ + fi'^dC) = 0}, which is orthogonal to spacelike hypersurfaces 



and has orbital angular velocity given by 
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This is the long-known Lense-Thirring precession [15| angular velocity arising due to the 
"dragging of inertial frame" effect of a stationary axisymmetric spacetime. Indeed, it is 
straightforward to demonstrate that this orthonormal tetrad observer can be identified with a 
ZAMO carrying zero intrinsic angular momentum. To this end, recall that when a spacetime 
metric possesses a rotational (azimuthal) isometry, there exists an associated rotational 
Killing field ip'^ = [d/dcf)^ = 5'^ such that its inner product with the tangent (velocity) 
vector = dx^ /dr (with r denoting the particle's proper time) of the geodesic of a test 
particle is constant along the geodesic, i.e., 

dt d(f) 

= 9,.J^ + 9,,J^- (10) 

Now, particularly when a local, stationary observer, which here is taken to be a test particle, 
carries zero angular momentum, L = 0, the angular velocity becomes 

_d(f) _ (dcp/dr) _ gt^ _ ^ 

dt {dt/dr) ■ ^ ' ^ ^ 



which confirms the identification of the local observer at rest in this orthonormal tetrad 
frame given above ZAMO. 

III. SOLUTION-GENERATING METHOD BY WALD 

1. Wald Field 

From the general properties of Killing fields Q| , a Killing vector in a vacuum spacetime 
generates a solution of Maxwell's equations in the background of that vacuum spacetime, 
long ago, Wald [3] constructed a stationary axisymmetric solution of Maxwell's equations 
in Kerr black hole spacetime. To be a little more concrete, Wald's construction is based on 
the following two theorems : 

(A) The axial Killing vector ip'^ = {d/dcj))^ generates a stationary axisymmetric test electro- 
magnetic field which asymptotically approaches a uniform magnetic field, has no magnetic 
monopole moment and has charge = AJ: 
= dip 



where "d" denotes the exterior derivative and J is the angular momentum of a Kerr black 
hole. 

(B) The time translational Killing vector = [d/dtY generates a stationary axisymmetric 
test electromagnetic field which vanishes asymptotically, has no magnetic monopole moment 
and has charge = —2M: 

where M is the mass of the Kerr hole. 

Apparently, this solution generating method of Wald can be applied in a straightforward 
manner to the construction of a stationary axisymmetric solution of Maxwell's equations in 
a spacetime surrounding a slowly-rotating neutron star represented by the Hartle-Thorne 
metric given earlier. Thus, we now look for the solution of the electromagnetic test field 
F that occurs when a stationary axisymmetric slowly-rotating neutron star of mass M and 
angular momentum J (represented by the Hartle-Thorne metric) is placed in an originally 
uniform magnetic field of strength Bq aligned along the symmetry axis of the neutron star. 
For now, we consider the case when the electric charge is absent. Obviously, then the solution 
can be readily written down, by referring to theorems (A) and (B) above, as follows : 
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Then, using F — \F^j,dx^ A dx"" and '0 = i^ndx'^, ^ — ^i^dx", with 
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and for the Hartle-Thorne metric given in the previous subsection, we can write the solution 
above in a concrete form as 
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Thus, this can be thought of as the Wald-type field for the case of a rotating neutron star. 
Of course, caution needs to be exercised, the point being that this Maxwell field is correct 
only for a slowly-rotating neutron star, but not for a neutron star with arbitrary rotation. 
2. Wald Charge 
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We now turn to the issue of charge accretion onto a slowly-rotating neutron star immersed 
in a magnetic field surrounded by an ionized interstellar medium ("plasma"). We shall 
essentially follow the argument again given by Wald j^, and to do so, we first need to know 
the physical components of the electric and the magnetic fields. This can be achieved by 
projecting the Maxwell field tensor given above in Eq. (14) onto the ZAMO tetrad frame. 
Thus, using the dual to the ZAMO tetrad, ca = (eo = e(t),ei = e(.r),e2 = e(e),e3 = e(0)), 
given earlier in Eq. 8), we can now read off the ZAMO tetrad components of F^i/ as 
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where (drg^), (drgu), (dggt^), (dggu), {drg^4,), and (dgg^^) are as given above in Eqs. (15)- 
(20). 

Here, consider particularly the radial component of the electric field (as observed by a 
local observer in this ZAMO tetrad frame) which, along the symmetry axis {6 = 0, vr) of 
the slowly-rotating neutron star, becomes 



Ef = El = F-] 
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Note that from this expression for the radial component of the electric field, it is not possible 
to determine the sign of charge accreted in terms of the relative direction between the angular 
momentum of the neutron star and the external magnetic field. This point is in contrast to 
what happens ,n the case of the Ke„. b.ae. hCeSf stud.ed p.eviouslv. To be a httle .o.e 
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concrete, it was realized that the radial component of the electric field for the case of Kerr 
black hole is radially inward/ outward if the hole's axis of rotation and the external magnetic 
field are parallel/ antiparallel. Put differently, this implies that if the spin of the hole and 
the magnetic field are parallel, then positively charged particles on the symmetry axis of the 
hole will be pulled into the hole while if the spin of the hole and the magnetic field are 
antiparallel, negatively charged particles on the symmetry axis of the hole will be pulled into 
the hole. In this manner, a rotating black hole will "selectively" accrete charged particles 
until it builds up an "equilibrium" net charge. For the present case of a rotating neutron 
star, however, an argument of this sort cannot be presented due essentially to the added 
complexity of the Hartle-Thorne metric compared to the Kerr metric. Nevertheless, one can 
still ask and answer the next natural question ; how then can the equilibrium net charge 
be determined ? To answer this question, we resort to the "injection energy" argument 
originally proposed by Carter for the Kerr black hole case. Again, it can be applied 
equally well to the present case of a slowly-rotating neutron star. 

Recall first that the energy of a charged particle in a stationary spacetime with time 
translational isometry generated by the Killing field = {d/dt)^ in the presence of a 
stationary electromagnetic field is given by 

e = -PaC = -gafSP'^f 

with = mu^ — eA^ being the 4-momentum of a charged particle with mass and charge ?fi 
and e, respectively. Now, if we lower the charged particle down the symmetry axis into the 
rotating neutron star, the change in electrostatic energy of the particle will be 

^ final ^initial ^^aC \horizon C^q,^ \oo 

Now, if < 0, it will be energetically favorable for the neutron star to accrete particles with 
this charge whereas if 6e > 0, it will accrete particles with opposite charge. In either case, 
the rotating neutron star will selectively accrete charge until is changed sufficiently that 
the electrostatic "injection energy" 6e is reduced to zero. We are, then, ready to determine, 
by this injection energy argument due to Carter, what the equilibrium net charge accreted 
onto the neutron star would be. 

In the discussion of the Wald-type field given above, we only restricted ourselves to 
the case of solutions to Maxwell equation in the background of an uncharged stationary 
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axisymmetric neutron star spacetime, and it was given by Eq. (12). Now, we need the 
solution when the stationary axisymmetric neutron star is shghtly charged via the charge 
accretion process described above. Then, according to theorem (B) in the discussion of 
the Wald field given earlier, there can be, at most, one more perturbation of a stationary 
axisymmetric vacuum (neutron star) spacetime which corresponds to adding an electric 
charge Q to the neutron star, which is nothing more than a linearly superposition of the 
solution {-Q/2M)F^ = {-Q/2M)d^ to the solution given in Eq. (12) to get 

F = -BoM + —di\ - —d^, (26) 
2 ^ M ^ 2M ^ ^ 

which, in terms of the gauge potental, amounts to 

A, = iMi^, + ^e.) - ^e.. (27) 

Then, the electrostatic injection energy can be computed as 
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where we used ^^^^ = —AR\r=ro,9=o,n (at r = ro, 6* = 0, vr), ^^^'^ = — 1 (at r — > oo, 6* = 0, vr), 
and ipfj,^'^ = (on the symmetry axis 6 = 0,tt) and tq denotes the radius of the neutron star. 
As was discussed carefully in Section II, vq is located outside of the gravitational radius of 
the rotating neutron star, which amounts, in our description, to the Killing horizon radius of 
the Hartle-Thorne metric at which A = 0. Thus, one may conclude that a rotating neutron 
star in a uniform magnetic field will accrete charge until the gauge potential evolves to a 
value at which 6e = 0, yielding an equilibrium net charge of Q = 2BqJ . Thus, this amount 
of charge may be called the Wald-type charge, and it is particularly interesting to note that 
it turns out to be the same as the "Wald charge" (the equilibrium net charge) for the case 

nn 

of a Kerr black hole p|, |9|] . This result appears to imply that the amount of equilibrium net 
charge Q = 2BqJ might be a generic value once a general relativistic rotating object gets 
charged via accretion. 

There, however, does exist a distinction between the case of a Kerr black hole and that 
of a rotating neutron star. That is, from Eq. (28), the rotating neutron star will eventually 
accrete an equilibrium net charge oi Q = 2BqJ only if [1 — AR\r=ro,e=o,Tr] 7^ 0. Certainly, 
this feature was absent in the case of a Kerr black hole and can be attributed to the different 
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structure of the Hartle-Thorne metric describing the region surrounding the rotating neu- 
tron star. Therefore, particularly if the values of mass (M), angular momentum (J), mass 
quadrupole moment (Q), and radius along the symmetry axis (ro) are such that they satisfy 

AR\r=ro,e=0,Tr (29) 



2M 2J2^ 

— + ^ 



1 Mr^ V ro/ 8 M ' 
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then it is "energetically unlikely" that a slowly-rotating neutron star would accrete any 
charge, and we find a peculiar feature which has no parallel in the Kerr black hole case. It, 
however, should be taken with some caution (i.e., should not be taken as being conclusive) 
because the injection energy argument due to Carter that leads to this conclusion is, rig- 
orously speaking, restricted to the situation along the symmetry axis 6* = 0, vr, and at this 
point it is not so obvious that the same would be true off the symmetry axis. 

Next, recall that we announced from the beginning that we shall consider the case when 
the charge accreted on the neutron star is small enough not to distort the background 
Hartle-Thorne geometry. Now, we provide the rationale that this is, indeed, what can 
actually happen. To do so, we first assume that just as in the case of a Kerr black hole 
(a particular slow- rotation limit of which is the Hartle-Thorne geometry), in the present 
Hartle-Thorne spacetime, J < M^. Then, since the typical value of the charge accreted on 
the rotating neutron star is Q = 2i?oJ, as just described, its charge-to-mass ratio has an 
upper bound 

i--«(i)^-«--(io^)(a-' 

where in the last equality we have converted and M from geometrized units to solar-mass 
n 

units and gauss 9]. Thus, for the typical case of a radio pulsar with mass M ~ Mq in a 
surrounding magnetic field of strength ~ 10^^(6"), QjM ^ 10^^ << 1. Thus, as we can 
see in this radio pulsar case, the charge-to-mass ratio of the associated rotating neutron star 
is small enough not to disturb the geometry itself. Thus, we can safely employ the solution- 
generating method suggested by Wald to construct the solution to Maxwell's equations when 
some amount of charge is around, to which we now turn. 

IV. STATIONARY AXISYMMETRIC MAXWELL FIELD AROUND A 
"SLIGHTLY CHARGED" AND SLOWLY-ROTATING NEUTRON STAR 
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As we discussed in the previous subsection, the stationary axisymmetric solution to 
Maxwell's equation in the background of a rotating neutron star with charge Q accreted 
in an originally uniform magnetic field can be constructed as 



(z2M)^d + 



-2M' 



Q 



2BnJ 



(30) 



Again for the Hartle-Thorne metric for a slowly-rotating neutron star, which was given in 
the previous subsection, the solution above can be written as 



(dog,,) + ^ (i - ^ 



{drgtt) 
(degtt) 
(drg^pt) 
I (deg^t) 



{dr A dt) 
{d9 A dt) 
{dr A d(f)) 
(d^ A#)], 



(31) 



and here, (drgtcp), {drgtt), {degt4>), (degtt), {drg^)^ {degn) ^ gi^^n earlier in the 
previous section. 

Obviously, in order to have some insight into the nature of this solution to Maxwell's 
equations, one may wish to obtain the physical components of the electric field and the 
magnetic induction. This can only be achieved by projecting the Maxwell field tensor above 
onto an appropriate tetrad frame like the ZAMO tetrad which, considering the Hartle- 
Thorne metric, can be well-defined as mentioned earlier. The ZAMO tetrad is well-known 
and is widely employed in various anayses in the literature. ZAMO is a fiducial observer 
following a timelike geodesic orthogonal to spacehke hypersurfaces. This implies that its 
4- velocity is just the timelike ZAMO tetrad u*^ = eg given earher in Eq. (8). Thus, if the 
Maxwell field tensor's components are projected onto the ZAMO tetrad frame, the physical 
components of the electric and the magnetic fields can be read off as Ei = Fio = F^j^(efeo) 
and Bi = eijkF^^ /2 = eyfcF'^^(e^e^)/2, respectively. Namely, from the dual to the ZAMO 
tetrad, ca = (eo = e(t), ei = e(r), 62 = e{e), 63 = e(^)), given earlier in Eq. (8), and 



'10 



Ax 1/2 
'S 



{AR)-'/'^Bo l^idrgt^) + 



2J 
M 



1 - 



Q 



2BoJ 



{drgtt) 
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20 



.^1/2 





+(Ai?)- 


1 

r3 2 


-^30 




= 0, 


Fl2 


— F p^p" 


= 0, 


Fi3 


— F P^p" 








1 



BoUdeg 



Q 



2J ( 



2BqJ 

Q \ 



25oJ, 



^ {degtt) 
(degc^t) \ 



(32) 
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F p^^p^ 



sin e 2 



^0 \ {drgn) + 



2J 

M 



Q 



1 



5, 







1 

M V 



25oJ 



{drg4,t 



r r/lV2 sin 2 - j^^-"-'^'^' ' m V" 2BoJ ^ 
V. MAGNETIC FLUX THROUGH A NEUTRON STAR 

With an asymptotically uniform stationary axisymmetric magnetic field aligned with the 
spin axis of a "slightly charged" slowly-rotating neutron star, as given above, we now would 
like to compute the flux of the magnetic field across one half of the surface of the neutron star, 
which is assumed to be of exact spherical geometry with radius tq. The physical motivation 
behind this study is to have some insight into the question of how much of the magnetic flux 
can actually penetrate the surface of a rotating neutron star - at least in idealized situations. 
A study of this sort might seem irrelevant since for a magnetized rotating compact object, 
all the flux would pass through its surface. Particularly for pulsars (which are believed 
to be rotating neutron stars), carrying intrinsic magnetic dipole moments, all the dipole 
magnetic field lines would be anchored on their surfaces. This naive picture, however, is 
based on electrodynamics in a flat spacetime and completely ignores the effects of highly 
self- gravitating rotating neutron stars. As we demonstrated earlier, in a region close to the 
surface of a rotating neutron star, Wald's solution exhibits possible higher-mult ipole fields, 
as well as a dipole field. Moreover, there is an interesting lesson we learned in the case of 
a Kerr black hole immersed in an asymptotically uniform magnetic field, and that is the 
observation that if the hole is maximally-rotating, the magnetic fields on the horizon get 
entirely expelled . Thus, one might wonder if the same can actually happen in the case of 
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rotating neutron stars. Indeed, these are the issues we would hke to address in this section 
in a careful manner. 

We now begin by considering two vectors lying on the surface of the neutron star. These 
vectors are given by 



dx'^ = (0, 0, de, 0), dx^ = (0, 0, 0, d(f)). 
Then, in terms of a 2nd-rank tensor constructed from these two vectors P|, 



rfo""^ = -{dx"dx2 — dxidx\ 



2)1 



one now can define the invariant surface element of the neutron star as 



(33) 



(34) 



ds = (2c/a„^(j"^)i/2 



1/2 



(35) 



Next, since the tensor dcr"^ is associated with the invariant surface element of any (not 
necessarily closed) 2-surface, the flux of the electric field and the magnetic field across any 
2-surfaces can be given, respectively, by 



(36) 



In particular, the flux of the magnetic field across one half of the surface of the neutron star 
is 



/■27r 


I-k/2 


/ 


h / dQFe^ \ 




Jo 



(37) 



Then, using Fq^ evaluated on the surface of the neutron star 



(38) 
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and the results of integration, Jq^'^ dO sin 6 cos = 1/2 and/o"/^d0(sin0 cos^6' — sin 0COS0) = 
0, we finally arrive at 



MroH ro y 8 M3 \ [r2(l - 2M/ro)]i/2'^2^M ' ^'^M 

This is the magnetic flux through one half of the surface of a slowly-rotating neutron star 
with accretion electric charge Q. Some discussions concerning the nature of this magnetic 
flux are now in order : 



(i) In the absence of accreted charge, i.e., (5 = 0, the magnetic flux above appears to drop 
to zero (i.e., the magnetic field gets entirely expelled) when the neutron star is maximally- 
rotating, which, based on the expression for the magnetic flux across the surface of the 
neutron star given in Eq. (39), amounts to J £„^ = Mtq/A. This curious feature, indeed, 
has a counterpart in the Kerr black hole case jlq], and that counterpart has been known 
for some time. For the present case of a neutron star, however, we are not fully qualified to 
make any definite statement on this point because the Hartle-Thorne metric employed rep- 
resents the region surrounding a slowly-rotating neutron star. Thus, the expression for the 
magnetic flux given above may be invalid for the case of rapidly-rotating neutron star, but 
that is beyond the scope of the present work as the associated metric is not available for now. 

(ii) For a nonvanishing accreted charge having a value in the range < Q < ^BqJ, it is 
interesting to note that 

$b(Q 7^ 0) > $b(Q = 0). (40) 

Actually, this is one of the points of central importance we would like to make in the present 
work. The physical interpretation of this characteristic can be briefly stated as follows : 
Once the neutron star accretes charge, toroidal currents occur as a result of rotation of the 
star. These toroidal currents, in turn, generate new magnetic fields which would be additive 
to the existing ones, increasing the total magnetic flux through the surface of the rotating 
neutron star. As in Ref. 9, essentially the same phenomenon takes place in the case of a 
slightly charged Kerr black hole. There, the interpretation was even clearer ; when the spin 
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of the hole and the asymtotically uniform magnetic field are parallel (antiparallel), the hole 
selectively accretes positive (negative) charge according to the injection energy argument 
proposed by Carter, and the charge, in turn, generates magnetic fields which are clearly 
additive to the existing ones. Also we expect that essentially the same will be true for the 
present rotating neutron star case although we could not determine the sign of the accreted 
charge in association with the relative direction between the spin of the neutron star and 
the direction of the magnetic field due to the highly complex structure of Hartle-Thorne 
metric, as we mentioned earlier. 

(iii) For the Wald charge Q = 2BqJ, the magnetic flux across the surface of the rotating 
neutron star appears to take on a rather special value : 

$B = Bonrl X (41) 
^MrlV^r-o) 8 \ [r2(l - 2M/ro)]V2'^2l^ ^) ^ ^) 

As mentioned in the introduction, it may be interesting to compare this analysis of the nature 
of the magnetic flux through a rotating neutron star with that of the Kerr hole studied in 
Ref. 9. That is, for all three cases (i), (ii), and (iii), the results discussed are generally 
consistent with those in the Kerr black hole case although the quantitative details are rather 
different between the two cases, the neutron star and the black hole, as the spacetime metrics 
representing the two are distinct. 

VI. CONCLUDING REMARKS 

In the present work, based on the solution-generating method given by Wald, the mag- 
netic fields around both uncharged and (slightly) charged neutron stars have been obtained. 
Particularly, for the charged neutron star, it has been demonstrated, following again the 
argument by Wald, that the neutron star will gradually accrete charge until it reaches an 
equilibrium value Q = 2BqJ. In association with the "magnetic braking" model of Goldre- 
ich and Julian P] for the spin-down of magnetized rotating neutron stars, the magnetic flux 
through one half of the surface of the rotating neutron star was computed, as well. For a 
nonvanishing accretion charge having a value in the range < Q < 2BqJ , the total magnetic 
flux through the neutron star has been shown to be greater than that without the accreted 
charge. Indeed, the physical interpretation of this characteristic can be briefly stated as 
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follows : When the spin of the neutron star and the direction of the asymtotically uniform 
magnetic field are parallel (antiparallel), the neutron star will selectively accrete positive 
(negative) charges following the injection energy argument proposed by Carter, and those 
charges, in turn, generate magnetic fields that add to the existing ones. 

Although our discussion thus far is equally relevant to all species of slowly-rotating rela- 
tivistic stars, indeed, we particularly had rotating neutron stars in mind. As such, it would 
be natural to attempt to make contact with the case of real pulsars as they are believed to 
be magnetized rotating neutron stars. Thus, below we consider how many of the theoretical 
results obtained in the present work can be carried over to a realistic, general relativistic 
description of the pulsar's magnetosphere. 

Radio pulsars are perhaps the oldest-known and the least energetic among all of the pulsar 
categories. The thoretical study of these radio pulsars can be traced back to the 1969 work of 
Goldreich and Julian ^. In their pioneering work, Goldreich and Julian argued that pulsars, 
which are thought to be the rotating magnetized neutron stars, must have a magnetosphere 
with a charge-separated plasma. They then demonstrated that an electric force, which is 
much stronger than the gravitational force, will be set up along the magnetic field and as 
a result, the surface charge layer cannot be in dynamical equilibrium. Then, there appear 
steady current flows along the magnetic field lines, which are taken to be uniformly rotating 
since they are firmly rooted in the crystalline crust of the pulsar surface. Although it was 
rather implicit in their original work, this model for the pulsar electrodynamics in terms of 
the force-free magnetosphere does, indeed, suggest that the luminosity of radio pulsars is due 
to the loss of their rotational energy (namely, the spin-down) and its subsequent conversion 
into charged particle emission, which eventually generates radiation in the far zone. 

Now, it is our hope to relate the results obtained in the present work to the physics of real 
pulsars ; thus, we begin with an inherent drawback of our treatment of the issues presented 
in this work. In the present work, we have mainly considered the case with a symmetric 
geometry in which the stationary axisymmetric magnetic field is precisely aligned with the 
neutron star's axis of rotation. Although this simplification/idealization was inevitable in 
order to make the system mathematically tractable, the magnetized rotating neutron star 
will fail to describe pulsars in the real world as the pulsars involve, by definition, oblique spin- 
magnetic field configurations. (In the case involving a rotating black hole, however, such 
an oblique configuration cannot be maintained because any possible initial misalignment 
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between the spin and the direction of the magnetic field will be removed in an efficient manner 
either by a viscous torque when the accretion disc is involved or by a magnetic alignment 
torque j20|.) Indeed, almost all of the theoretical studies of pulsar electrodynamics assume 
this idealized symmetric geometry for concrete, analytical calculations, and concerning this 
limitation, our work makes no progress. Thus we are still far away from a satisfactory 
theoretical description of the real physics that takes place in the pulsar's magnetosphere. 
Nevertheless, there are quite a few features represented by the theoretical results in the 
present work that could be relevant to the real pulsar situation. To make this point indeed 
look convincing, we first comment on our motivation for choosing the Wald algorithm to 
obtain the solutions of Maxwell's equations in the background spacetime of neutron stars. 
Once again, the philosophy employed in the present work is to take the electromagnetic 
fields surrounding a rotating neutron star as a theoretical model to represent the pulsar 
magnetosphere (except possibly for the plasma content). Therefore, as a theoretical model, 
our approach is, indeed, not without some inherent limitations. We here remind the reader 
that earlier in the introduction, we mentioned that the results presented here for the case 
of a rotating neutron star can be compared in parallel with those that have been presented 
in our earlier work [9] for the case of a rotating Kerr black hole because we employed 
in this work the same algorithm of Wald to determine the structure of (electro)magnetic 
field around the rotating neutron star. Now, considering this point, notice that in Wald's 
construction of a stationary axisymmetric solutions of Maxwell's equations in a Kerr black 
hole spacetime or in a Hartle-Thorne spacetime for a rotating neutron star, we assumed 
on technical grounds from the outset that the Kerr black hole or the rotating neutron star 
was placed in an originally uniform magnetic field. Indeed, this setup could be relevant 
for objects like pulsars which produce very strong external dipole magnetic fields. After 
all, this asymptotically uniform magnetic field configuration is, at least, not inconsistent 
with the open dipole magnetic field lines in the far zone in the "aligned rotator" model 
of Goldreich and Julian p]. Then, one looks for a surrounding (electro) magnetic field, 
which gets distorted by the presence of a Kerr black hole or a rotating neutron star. In 
a rigorous sense, therefore, this boundary condition in Wald's construction may not seem 
very relevant for real pulsar situations in which pulsars are believed to carry dipole moments 
and, hence, produce dipole magnetic fields. In this regard, our approach based on Wald's 
construction method can be viewed as one seeking systematic gravitational corrections to an 
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intrinsic dipole magnetic field. Indeed, as we stated earlier, Wald's solution has no monopole 
moment. In the region close to a neutron star, however, it may well have a non-trivial dipole 
and possibly higher mutipole fields. Also, we beheve that this behavior of Wald's solution is 
a new contribution to the existing information on the structure of a pulsar magnetosphere 
given in the literature. 

Lastly, one might be worried about the validity of the Hartle-Thorne metric for the region 
surrounding the slowly-rotating neutron stars employed in this work to describe the mag- 
netosphere of pulsars, which seem to be rapidly-rotating and to have typically millisecond 
pulsation periods. Thus, in the following, we shall defend this point in a careful manner. 
Here, "slowly-rotating" means that the neutron star rotates relatively slowly compared to 
an equal-mass Kerr black hole which can rotate arbitrarily rapidly up to its maximal rota- 
tion of J = M^. Thus, this does not necessarily mean that the Hartle-Thorne metric for 
slowly-rotating neutron stars cannot properly describe millisecond pulsars. To see this, note 
that according to the Hartle-Thorne metric, the angular speed of a rotating neutron star 
is given by the Lense-Thirring precession angular velocity in Eq. (9) at the surface of the 
neutron star, which, restoring the fundamental constants to get back to Gaussian units, is 



As we mentioned earlier, one of the obvious differences between the black hole case and the 
neutron star case is the fact that the black hole is characterized by its event horizon while 
the neutron star has a hard surface. As such, in terms of the spacetime metric generated by 
each of them, just as the Lense-Thirring precession angular velocity (due to frame-dragging) 
at the horizon represents the black hole's angular velocity the Lense-Thirring precession 
angular velocity at the location of the neutron star's surface should give the angular velocity 
of the rotating neutron star. 

Thus, the Hartle-Thorne metric gives the angular speed of a rotating neutron star having 
the parameters of a typical radio pulsar, M ~ 2 x 10^^ (5^) and Tq ~ 10^ (cm), as a; = 
2a{M'^ /r^){G/c){G/c'^) ~ 10^(l/sec) which, in turn, yields a rotation period of r = 27i/uj ~ 



10-2(sec). Here we used (G/c^) = 0.7425 x 10-^^{cm/g) and (G/c) = 2.226 x 10-^^{cmyg- 



sec). Indeed, this is impressively comparable to the observed pulsation periods of radio 
pulsars, T ~ 10~^ — l(sec), as we discussed earlier. As a result, we expect the Hartle-Thorne 
metric to be well-qualified to describe the geometries of millisecond pulsars. 




(0 < a < 1). 



(42) 
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To summarize, some aspects of the results of the theoretical study of a rotating neutron 
star's electrodynamics carried out in this work seem irrelevant while some appear to be quite 
relevant to real pulsar situations. Namely, although the predictions of the present theoretical 

study of a rotating neutron star's electrodynamics are not fully satisfying, it is, nevertheless, 
our hope that at least here we have taken one step closer to a systematic general relativistic 
study of the electrodynamics in a region close to a rotating neutron star in association with 
their pulsar interpretation. 
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